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A basic prerequisite for in vivo X-ray imaging of the lung is the exact 
determination of radiation dose. Achieving resolutions of the order of 
micrometres may become particularly challenging owing to increased dose, 
which in the worst case can be lethal for the imaged animal model. A framework 
for linking image quality to radiation dose in order to optimize experimental 
parameters with respect to dose reduction is presented. The approach may find 
application for current and future in vivo studies to facilitate proper experiment 
planning and radiation risk assessment on the one hand and exploit imaging 
capabilities on the other. 



1. Introduction 

The development of high-speed time-resolved tomographic 
microscopy is of great interest for various three-dimensional 
in vivo studies. One important apphcation is the study of lung 
dynamics, in particular lung inflation and deflation issues 
during physiological as well as mechanical ventilation, which is 
required, for example, after a premature birth or during a 
general anesthetic. Two hypotheses on the structural altera- 
tions in the gas-exchange area during breathing are still under 
debate: a heterogeneous distention pattern (Mertens et al, 
2009) of different lung areas and a homogeneous cycHc 
opening-and-coUapse (Albert et al. , 2009) of all alveoli. 

In the recent past, lung imaging with small animal models 
has become an estabhshed technique at synchrotrons (Yagi et 
al, 1999; Bayat et al, 2001; Lewis et al, 2005), making rapid 
data acquisition and high-resolution X-ray imaging possible 
(Flannery et al., 1987). Kitchen et al. (2004) performed phase- 
contrast X-ray imaging of mice and rabbits in vivo in two- 
dimensions and developed a method for visualizing lung liquid 
clearance at birth (Kitchen et al., 2005; Hooper et al, 2009). 
Bayat et al. (2001) used xenon as a contrast agent and were 
able to measure and visualize gas distribution within the lungs 
in three dimensions, enabling further functional studies 
(Suhonen et al, 2008; Bayat et al, 2008, 2009). More recently, 
synchrotron-based tomographic microscopy has been utilized 
to study the distribution of gas flow throughout the airway tree 
in connection with altered lung motion as an indicator for 
regional lung disease (Fouras et al, 2012; Dubsky et al, 2012). 

So far, all studies either were performed in two dimensions 
only or suffered from low temporal and spatial resolution; 
high-resolution lung images in three dimensions, on the other 



hand, were successfully obtained only for static samples 
(Schittny et al, 2008; Yong et al, 2009; Haberthur et al, 2010; 
Mokso et al, 2011; Zhang et al, 2011). Thus, in vivo X-ray 
tomography with micrometre spatial and sub-second temporal 
resolution remains a challenge and many open questions in 
the study of lung physiology remain unanswered. 

We describe our approach to image formation of biologi- 
cally relevant features in the lung, aiming at optimal image 
quality in terms of contrast, spatial and temporal resolution, 
and deposited radiation dose. In particular, we show current 
limitations towards in vivo imaging at the micrometre scale. 



2. Materials and methods 
2.1. Image acquisition 

The experiment was carried out at the X02DA TOMCAT 
beamline of the Swiss Light Source (SLS) at the Paul Scherrer 
Institute (Villigen, Switzerland). The experimental setup is 
depicted in Fig. 1: the X-ray beam, produced by a 2.9 T 
bending magnet on a 2.4 GeV storage ring (with ring current 
/ = 400 mA, top-up mode), is monochromated with a double- 
multilayer monochromator and tuned to 21 keV. 

A sample-to-source distance of 25 m is used for producing 
an X-ray beam with appropriate spatial coherence properties. 
We used a high-speed CMOS detector (PCO.Dimax) coupled 
to visible-hght optics with 100 and 20 |im-thick scintillators for 
medium and high spatial resolutions, respectively. The samples 
were probed with two different optics, yielding effective pixel 
sizes of 2.9 and 1.1 jim, respectively. For these two optics the 
field of view was adjusted with horizontal and vertical slits. 
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located just before the sample and producing beam sizes of 
5.8 X 2.7 mm and 2.2 x 2.2 mm, respectively. The sample-to- 
detector distance z was varied in the range 24-300 mm to 
allow for the variation of the Fresnel interference pattern for 
image quality optimization purposes. Raw images were 
acquired with exposure times ranging from 2 to 13 ms and 901 
tomographic projections. 

The measurements were performed ex vivo on mice aged 
37 d (/t = 3 / Balb-C, central animal facility of the University of 
Bern) that were sacrificed before the experiment. The mice 
were killed with an overdose of a combination of Acepro- 
mazine, Xylazin and Ketamin. They were then placed in an 
upright position into a 2.5 cm-diameter Falcon tube. All parts 
of the animal experiments were approved and supervised by 
the Swiss Agency for the Environment, Forest and Landscape, 
and the Veterinary Service of the Canton of Bern. 



2.2. Post-processing 

The aforementioned setup facilitates propagation-based 
phase contrast, represented by interference fringes on the 
tissue interfaces in the recorded digital projection images. 
Prior to further analysis, each projection was corrected with 
the respective dark and flat-field image. In a second step, 
single-image phase and intensity extraction were applied to all 
projections (Paganin et al, 2002). In a third step, computed 
tomography (CT) reconstruction was conducted with the 
gridrec algorithm (Marone & Stampanoni, 2012), enabling fast 



reconstructions of large data sets. Finally, the CT reconstruc- 
tions obtained from phase-retrieved images and absorption 
images were fused in the Fourier domain in order to correct 
for the high-pass characteristics of the Paganin algorithm 
(Irvine et al., 2013). The complete post-processing flowchart is 
depicted in Fig. 2. 

2.3. Image analysis 

The reconstructed tomographic slices were examined in 
view of contrast-to-noise ratio (CNR) and resolution. For 
CNR the following formula was used: 



CNR = 2 
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Figure 1 

Experimental setup at the X02DA TOMCAT beamline. 
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Figure 2 

Flowchart of the post-processing pipeline. After image acquisition, the CT-reconstructed volumes 
obtained from absorption and phase-retrieved images are fused in the Fourier domain. Afterwards, 
further analysis (segmentation, visualization etc.) can be apphed. 



where Si represent mean pixel values and cr, standard devia- 
tions of a manually defined object and background region, 
respectively. One has to consider, however, that the object (i.e. 
lung tissue) consists of blood vessels and erythrocytes which 
have different densities and may lead to distorted contrast 
levels. Under the present imaging scheme, this effect is 
negligible. 

The resolution was determined by a criterion based upon 
Fourier analysis (Modregger et al, 2007). Considering a test 
image the resolution was obtained by taking a line profile and 
calculating its power spectral density (PSD). The PSD 
converges towards a value that can be defined as the noise 
baseline. Taking the value of the PSD at the noise baseline 
twice and matching it to the respective 
spatial frequency yields the resolution. 
This can be made more robust and fully 
automated if one operates on the mean 
PSD, obtained from many line profiles. 
The only necessary condition for this 
method requires the highest spatial 
frequency of the image to be lower than 
the Nyquist frequency of the line 
profile. Since the smallest features of 
the lung are at least of the order of tens 
of pixels, this condition was fulfilled for 
both optics of our setup. 

Additionally, the results were cross- 
checked with another method by which 
the resolution is derived from a line 
profile taken along an edge in the image. 
The line profile is then fitted with an 
error function by means of a least- 
squares fit and the resolution is deter- 
mined from the slope of the function. 

In Fig. 3 the test image with the 
defined object and background regions, 
the line profile, the line profile along the 
edge, and the calculated PSD are 
shown. The image in Fig. 3(a) was 
constructed from a region near the 
center of the original tomographic slice 
in order to reduce the radial depen- 
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dency of the resolution originating from the tomographic 
reconstruction. The resolution Xres is calculated by 



•^res /^sizeC-^n/^res)' 



(2) 



where p^i^^ is the pixel size of the detector, x„ the number of 
pixels for the taken line profile and A:res the spatial frequency 
obtained from the resolution criterion in Fig. 3 (J). 

2.4. Dose calculations 

The X-ray radiation dose was assessed in two steps by 
theoretical and experimental means. First, the theoretical flux 
F was calculated with the TOMCAT beamline parameters 
(Stampanoni et al, 2006) and the following formula (Kim, 
1995): 
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Figure 3 

{a) A tomographic slice of a lung for the 2.9 )j,m-pixel-size optics. Regions 
for calculating CNR are defined by the rectangles, the red dashed line 
indicates the input for the line-profile plot in {b) and the solid magenta 
line indicates the edge line plot in (c). {d) Mean PSD from many line plots 
with the calculated noise baseline and the respective resolution criterion. 




Figure 4 

Dose and flux measuring setup, {a) Water column with reference 
ionization chamber and a silicon pin diode (placed behind the ionization 
chamber, {b) Frozen mouse sample. 
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where cOc is the critical frequency of the bending-magnet 
radiation, co is the emitted photon frequency and 7^5/3 repre- 
sents the modified Bessel function of the second kind. In 
equation (3) the flux is given in practical units [photons 
s~^mrad~^ (0.1%, bandwidth)"^]. Thus it has to be further 
corrected in order to include bandwidth and reflectivity of the 
monchromator as well as ah optical elements that further 
reduce the theoretical flux. In a second step, the absorbed dose 
(Asadchikov et al., 2010) of a 2.5 cm water column was 
calculated from NIST mass attenuation coefficients (http:// 
physics.nist.gov/PhysRefData/XrayMassCoef/ComTab/water. 
html). 

The X-ray flux was determined experimentally with a cali- 
brated silicon pin diode (Owen et al, 2009) and a frozen 
mouse sample in connection with a 2.5 cm water column to 
verify that the modeling with water is appropriate. Images 
from the setup are depicted in Fig. 4. 

3. Results 

Lung images were reconstructed from tomographic projec- 
tions taken at five different propagation distances. Fig. 5 shows 
the calculated CNR and resolution for each distance and for 
different weighting factors for the absorption images in the 
fusion algorithm. The values in brackets denote resolutions 
obtained from the edge-fitting method. As expected, for the 
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Figure 5 

Contrast-to-noise ratio (CNR) as a function of sample-to-detector 
(propagation) distance and the weighting of the absorption image in 
the fusion algorithm. The values were calculated for the 2.9 )im-pixel-size 
optics and 11.7 s total scan time. The values in brackets denote 
resolutions obtained from the edge-fitting method. 
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transport-of-intensity-based phase-retrieval approach in the 
near-field region the optimal propagation distance is near the 
value of pliJ)^ (m), where /7size is the pixel size and A is the 
wavelength. We show that fine tuning this distance can 
significantly affect CNR. By applying the fusion algorithm 
(Irvine et aL, 2013), it is possible to increase the resolution 
further, but at a cost of CNR. This can be used for adjusting 
the resolution of the reconstructed image. 

It has been reported elsewhere (Lewis, 2004; Sera et al., 
2008) that phase-contrast images of lungs exhibit better 
contrast and especially less noise than pure absorption images. 
We addressed this issue by evaluating the degree of 
'segmentability' in a tomographic scan. In Fig. 6, two images 
and the respective binary images obtained by thresholding are 
plotted. The total scan time was 5.4 s and the CNR values for 
the absorption and phase-retrieved images were 1.8 and 15.0, 
respectively, which makes up a difference of approximately 
one order of magnitude. However, the segmentation from the 
absorption images shows small artefacts originating from the 
higher noise level. We further found that a CNR of ~2 is 
currently the lower limit for a successful segmentation and in 
some cases necessitates further processing (e.g. filtering). This 
limit was verified by visual inspection. 

Finally, CNR and the X-ray radiation dose were calculated 
as functions of the total scan time of a tomographic scan, as 
shown in Fig. 7. The values for CNR obtained from experi- 
mental data were fitted and extrapolated using the function 



f(t) = a[l - exp(-bt)l 



(5) 



where a, b represent arbitrary parameters for the nonlinear fit 
and t the total scan time. The function derives from the fact 
that CNR will be zero if the sample is not exposed to X-rays 
(t = 0). Thus its origin is at zero. The upper value of CNR is 
limited by the mean pixel values represented by the density 




Figure 6 

Comparison of tomographic slices before and after thresholding obtained 
from pure absorption-contrast (a) and phase-retrieved projection images 
(b). For both images the total scan time was 5.4 s. The small rectangles 
denote the regions for calculating CNR. 



values of the biological material and their standard deviations, 
which are inherent to the acquisition scheme. CNR will 
therefore converge to a saturation value in the tomographic 
reconstruction. 



4. Discussion 

The optimization of CNR with respect to the deposited 
radiation dose is crucial in low-dose experiments. In particular, 
we found that the application of the one-shot phase-retrieval 
algorithm (Paganin et al., 2002) increased CNR by almost a 
factor of ten compared with the pure absorption-contrast 
images. However the spatial resolution decreased, as shown in 
Fig. 5. It is also evident that tuning the propagation distance 
can affect CNR by a factor of two under the given conditions. 
These results indicate that the optimal propagation distance, 
represented by a trade-off between CNR and resolution, and 
fusion parameters are particularly important and should be 
considered before the experiment. Furthermore, we found 
that the resolutions obtained from the two different criteria 
are not consistent for larger propagation distances and the 
robustness of the edge-fitting method in general was bad. This 
can be explained as follows: by increasing the propagation 
distance the edge enhancement becomes more pronounced 
and the Fourier resolution criterion detects edge artefacts as 
'features', whereas the edge fitting gives a poorer result. Thus, 
the Fourier resolution criterion has to be handled with care for 
larger propagation distances. 

For the radiation dose we found that the 1.1 |im-pixel optics 
(numerical aperture: 0.4) yielded a lower dose than the 
2.9 |im-pixel optics (numerical aperture: 0.2) for the same 
CNR values. The main reason for this result is the fact that 
images for the two optics were not acquired under the same 
experimental conditions, i.e. optimized propagation distance 
etc. Apart from that, one also has to regard the efficiencies of 
the different optics. Finally, another reason is the fact that 
images from both optics were reconstructed from local 
tomographic projections. Since we calculated the total dose in 
the irradiated volume, the smaller field-of-view means lower 
energy deposition for the same volume. Obviously, the dose 
rate in the lungs will be higher for the 1.1 |im-pixel optics. 
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Figure 7 

CNR as a function of total scan time per tomographic scan for the two 
different optics in use. (a) Optics with 2.9 [im effective pixel size, (b) 
Optics with 1.1 \xm pixel size. 
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The model function from equation (5) for extrapolating 
CNR as a function of the total scan time still needs to be 
verified with the gridrec algorithm (Marone & Stampanoni, 
2012). However, taking the range of interest that we analyzed, 
a similar fit (e.g. linear) would only cause negligible errors. 
From the results in Figs. 6 and 7 we further hypothesize that 
under the given conditions the lowest achievable dose for 
obtaining the necessary CNR is 10 and 5 Gy, respectively. 

5. Conclusion 

We have developed a framework for optimizing experimental 
parameters in medium-dose experiments, which is particularly 
important for current efforts in developing microtomographic 
in vivo X-ray imaging to study lung physiology at the micro- 
metre scale. Thanks to the fast phase-retrieval and CT 
reconstruction algorithms, it is possible to apply the optimi- 
zation steps on-the-fly (or even in real time) as an initial part 
of a beamline experiment. 

We showed that the lowest achievable dose at the moment 
is in the range of 5-10 Gy per tomographic scan at a total scan 
time of approximately 0.5 s and approximate resolutions 
between 4-10 |im, producing images with an approximate 
CNR of 2. Our results indicate that in vivo tomography at the 
micrometre scale and sub-second temporal resolution should 
be feasible, but will necessitate further adjustments. 

We thank Gordan Mikuljan for his great support in setting 
up the experiment at the beamline. This study was supported 
by the National Competence Center in Biomedical Imaging 
(NCCBI), grant No. 1126.0076, and the Swiss National Science 
Foundation (SNF), grant No. CR23I2-135550. 
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